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ABSTRACT 

In  a  recent  paper  [ Phys .  Fluids  .42:1678-1684,  1990],  the  authors  compared  the  perfor¬ 
mance  of  a  variety  of  turbulence  models  including  the  K  —  e  model  and  the  second-order 
closure  model  derived  by  Yakhot  and  Orszag  based  on  Renormalization  Group  (RNG)  meth¬ 
ods.  The  performance  of  these  RNG  models  in  homogeneous  turbulent  shear  flow  was  found 
to  be  quite  poor,  apparently  due  to  the  value  of  the  constant  (7ei  in  the  modeled  dissipation 
rate  equation  which  was  substantially  lower  than  its  traditional  value.  However,  recently  a 
correction  has  been  made  in  the  RNG  based  calculation  of  Ce i-  It  is  shown  herein  that  with 
the  new  value  of  Ce\,  the  performance  of  the  RNG  K  —  e  model  is  substantially  improved. 
On  the  other  hand,  while  the  predictions  of  the  revised  RNG  second-order  closure  model  are 
better,  some  lingering  problems  still  remain  which  can  be  easily  remedied  by  the  addition  of 
higher-order  terms. 


’This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No.  NAS1-18605  while  the  first  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 
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A  comparative  study  of  the  performance  of  nine  independent  turbulence  models  in  ro¬ 
tating  homogeneous  shear  flow  was  recently  reported  by  Speziale  et  al.1  Two  of  the  models 
considered  consisted  of  the  K  —  e  model  and  second-order  closure  model  derived  by  Yakhot 
find  Orszag2  using  Renormalization  Group  (RNG)  methods.  It  was  rather  surprising  how 
poorly  the  RNG  models  performed  in  homogeneous  shear  flow  relative  to  the  older,  em¬ 
pirically  based  models  of  the  same  general  type.  The  origin  of  the  deficient  predictions  of 
the  RNG  models  appeared  to  be  largely  due  to  the  rather  low  value  of  the  constant  C, i  in 
the  modeled  dissipation  rate  equation;  the  RNG  value  of  Ct i  was  1.063  in  contrast  to  the 
more  traditional  value  of  Ct\  =  1.44.  However,  a  recent  re-examination  of  the  RNG  based 
calculation  of  Ct\  by  Yakhot  and  Smith3  has  led  to  a  correction  -  the  new  value  of  Cc\  is 
1.42.  Some  minor  changes  in  the  values  of  other  constants  in  the  RNG  K  —  e  model  were  also 
made.3  In  light  of  these  changes,  it  would  be  desirable  to  set  the  record  straight  in  regard 
to  what  these  Renormalization  Group  models  now  predict  for  homogeneous  shear  flow  -  a 
critical  test  case  used  to  evaluate  the  performance  of  models.  This  establishes  the  motivation 
for  the  present  paper. 

In  the  RNG  K  —  e  model,  the  Reynolds  stress  tensor  r,j  =  u'u'  (given  that  u't  is  the 
fluctuating  velocity  find  an  overbar  represents  an  ensemble  mean)  is  modeled  as  follows:2’3 


TH 


(1) 


where  K  =  \u'xu[  is  the  turbulent  kinetic  energy,  e  =  vdu'Jdxj  du[/dxj  is  the  turbulent  dis¬ 
sipation  rate,  Ui  is  the  mean  velocity,  and  is  a  dimensionless  constant  which  is  calculated 
to  be  0.085.  In  homogeneous  turbulence,  the  turbulent  kinetic  energy  is  a  solution  of  the 
transport  equation 


•  dui 

K  =  -t”6^-S 


(2) 


which  is  exact.  The  turbulent  dissipation  rate  is  obtained  from  the  RNG  derived  transport 
equation 
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where  Ct\  —  1-42  and  Ce2  =  1.68  according  to  the  recent  calculations  of  Yakhot  and  Smith.3 
These  new  values  constitute  a  correction  to  the  earlier  values  of  Ce i  =  1.063  and  =  1.72 
reported  by  Yakhot  and  Orszag.2  An  additional  production  term  was  also  uncovered  by 
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Yakhot  and  Smith3  which  they  were  unable  to  close.  However,  an  order  of  magnitude 
analysis3  indicated  that  this  term  is  small  unless  there  are  large  strain  rates  -  a  case  which 
will  not  be  considered  herein.  Hence,  we  will  neglect  this  additional  term  in  the  present 
study.  For  the  RNG  second-order  closure  model,  the  eddy  viscosity  model  (1)  is  replaced 
with  a  Reynolds  stress  transport  model  of  the  form1 


where  C\  and  Cj  are  constants  that  are  calculated  to  be  1.59  and  2/15,  respectively.  Some 
clarifications  are  needed  concerning  the  origin  of  this  model  which  has  not  been  published 
and  was  obtained  from  a  private  communication  with  V.  Yakhot.  We  have  come  to  learn 
that  this  was  not  intended  to  be  a  final  model,  but  rather  was  the  result  of  a  low-order 
calculation  of  the  pressure- strain  correlation  whose  purpose  was  to  merely  demonstrate  that 
the  Rotta  term  -  with  a  coefficient  C\  close  to  the  well  accepted  value  of  1.5  -  could  be 
formally  obtained  from  RNG.  Hence,  the  results  predicted  by  this  preliminary  model  should 
be  judged  accordingly. 

In  homogeneous  shear  flow,  an  initially  isotropic  turbulence  where 
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Tij  =  — A"o£,'yj  £  =  Eq  (5) 


at  time  t  =  0  is  subjected  to  a  constant  shear  rate  S  with  the  corresponding  mean  velocity 


gradient  tensor 


dui 

dxj 
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In  Figure  1,  the  time  evolution  of  the  turbulent  kinetic  energy  (where  K*  =  K/Kq  and 
t*  =  St)  predicted  by  the  new  RNG  K  —  e  model  is  compared  with  the  large-eddy  simulation 
of  Bardina  et  al .*  for  an  initial  condition  of  e0/SKo  =  0.296.  The  predictions  of  the  old 
version  of  the  RNG  K  —  s  model  (where  C ^  =  0.0837,  Ct\  =  1.063,  and  Gt 2  =  1.72)  as  well 
as  the  standard  K  —  £  model  (where  C =  0.09,  Ce\  =  1.44  and  Ce a  =  1.92)  are  also  shown 
in  Figure  1.  It  is  clear  from  these  results  that  the  revised  RNG  K  —  e  model  does  the  best 


overall  job  in  reproducing  the  growth  rate  of  the  numerical  experiment  on  homogeneous  shear 
flow.  Analytically,  it  can  be  shown  why  this  is  the  case.  From  a  straightforward  calculation, 
it  can  be  shown  that  the  turbulent  kinetic  energy  and  dissipation  rate  grow  exponentially  in 


homogeneous  shear  flow  as  follows:1’5 

K*  ~  exp(Af*),  e*  ~  exp(At*) 
where  the  dimensionless  growth  rate  A  is  given  by 


A  = 


C^Ct2  -  Cei)2 


-.1/2 


(Ctl  -  l)(Ce2  -  1) 


(7) 


Hence,  the  growth  rate  becomes  singular  when  Cei  =  1  —  a  state  of  affairs  that  explains 
why  the  old  version  of  the  RNG  K  —  e  model,  with  Ce i  =  1.063,  overpredicted  the  growth 
rate  of  the  turbulent  kinetic  energy  by  such  a  wide  margin.  The  new  version  of  the  RNG 
K  —  e  model  predicts  a  growth  rate  of 


A  =  0.142 


which  is  extremely  close  to  the  range  of  values  obtained  from  physical  and  numerical 
experiments.6,7  On  the  other  hand,  the  standard  K  —  e  model  predicts  the  somewhat  high 
value  of  A  =  0.226  which  explains  why  this  model  overpredicts  the  LES  data  for  K*  as  shown 
in  Figure  1 .  A  more  complete  set  of  the  equilibrium  values  predicted  by  these  different  ver¬ 
sions  of  the  K  —  £  model  will  be  provided  later. 

In  Figure  2,  the  time  evolution  of  the  turbulent  kinetic  energy  predicted  by  the  revised 
RNG  second-order  closure  model  is  compared  with  the  large-eddy  simulation  of  Bardina  et 
al*  as  well  as  with  the  predictions  of  the  earlier  version  of  the  model  and  the  Launder,  Reece, 
and  Rodi®  (LRR)  model.  The  new  version  of  the  RNG  model  does  yield  better  predictions 
than  the  older  version  of  the  model  since  the  previous  value  of  Cel  =  1.063  was  too  close 
to  C€\  =  1  which  constitutes  a  bifurcation  point  of  the  dissipation  rate  transport  equation 
as  shown  by  Speziale.9  However,  there  are  still  problems  with  the  model  which  gives  rise  to 
points  of  inflection  in  the  time  evolution  of  K*  -  a  feature  that  makes  it  inferior  ‘o  other 
second-order  closure  models  such  as  the  Launder,  Reece,  and  Rodi  model.  The  origin  of 
this  problem  appears  to  be  tied  to  the  modeling  of  the  pressure  strain  correlation.  In  the 
Launder,  Reece,  and  Rodi  model,  the  pressure  strain  correlation  n„  =  p'(du'Jdxj  +  du'j/d Xi) 
is  modeled  as  follows 


nu  —  —2Cicbtj  -f  2C2K$ij  +  C3K  (bikSjk  +  bjk^ik  ~  qbkiS uSijj 


+C<K(bikWjk  +  bjkWik) 


(8) 
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where 


(in  the  simplified  form  of  the  Launder,  Reece,  and  Rodi  model,  C\  =  1.8,  C2  =  0.4  and 
C3  =  C4  =  1.2).  This  model  satisfies  two  important  consistency  conditions:  (a)  the  constant 
C2  is  equal  to  2/5  -  a  result  that  follows  from  simple  symmetry  conditions  for  Iljj  as  well 
as  from  Rapid  Distortion  Theory,8  and  (b)  it  represents  a  formal  expansion  of  Ifr,-  to  0(b) 
in  the  anisotropy  tensor.  On  the  other  hand,  for  this  preliminary  RNG  second-order  closure 
model  we  have 

Uij  =  -2  Cxebij  +  2C2Kl>ij  (11) 

where  C\  =  1.59  and  C2  =  This  model  is  not  complete  to  0(b)  in  the  rapid  pressure- 
strain  term  and  violates  the  important  symmetry  constraint  of  C2  =  2/5.  The  fixed  points 
that  the  resulting  nonlinear  ODE’s  for  these  second-order  closure  models  give  rise  to  in 
homogeneous  shear  flow  are  of  the  focus  type.5  Significant  deviations  of  C2  from  a  value 
of  2/5  excites  the  imaginary  parts  associated  with  these  fixed  points,  thus  inducing  inertial 
oscillations  which  are  unphysical  for  the  case  of  pure  shear  flow. 

An  overview  of  the  performance  of  the  models  can  be  gleaned  from  Table  1  which  com¬ 
pares  the  predicted  equilibrium  values  with  the  most  recent  experimental  data  of  Tavoularis 
and  Karnik7  for  homogeneous  shear  flow  (this  data  constitutes  a  mean  over  the  stronger 
shear  rate  cases).  Here  (•)„,  denotes  the  equilibrium  value  obtained  in  the  limit  as  t  — *  oo. 
Several  observations  concerning  Table  1  are  noteworthy: 

(a)  The  revised  RNG  K  —  e  model  yields  substantially  better  results  than  the  old  version 
of  the  model  and  is,  on  balance,  better  than  the  standard  K  —  e  model.  This  appears  to 
explain  why  the  models  performed  as  they  did  in  Figure  1  relative  to  the  LES  results. 

(b)  The  only  deficiency  in  the  predictions  of  the  new  RNG  K—e  model  for  homogeneous  shear 
flow  are  in  the  values  of  the  normal  components  of  the  anisotropy  tensor  -  a  shortcoming  of 
any  model  based  on  an  isotropic  eddy  viscosity.  However,  the  RNG  based  anisotropic  eddy 
viscosity  model  of  Rubinstein  and  Barton10  -  which  predicts  (&n)oo  =  0.260  and  (&2a)oo  = 
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—0.196  for  the  normal  anisotropies  in  homogeneous  shear  flow  -  alleviates  this  deficiency  to 
a  large  extent. 

(c)  The  RNG  second-order  closure  model  does  perform  somewhat  better  with  the  new  value 
of  Ct i  (the  model  now  predicts  a  weak  exponential  time  growth  of  K*  whereas  the  old 
version  of  the  model  predicted  a  power  law  growth,  with  ( SK/e)oo  =  oo,  due  to  the  close 
proximity  of  Ce i  to  the  bifurcation  point  Ce i  =  1).  However,  this  preliminary  model  still 
performs  weakly  in  comparison  to  the  more  commonly  used  second-order  closures  such  as 
the  Launder,  Reece,  and  Rodi  model.  The  deficiency  in  this  model  is  traced  to  the  rapid 
part  of  the  pressure-strain  correlation  which  is  0(1)  instead  of  0(b)  in  the  anisotropy  tensor. 
In  fact,  the  deviation  of  C2  from  0.4  to  2/15  results  from  the  model  trying  to  compensate 
for  the  truncated  0(b)  terms11  (interestingly  enough,  if  C2  is  set  to  0.4  in  Eq.  (11),  the 
predictions  of  the  model  deteriorate  substantially).  Hence,  we  have  little  doubt  that  if  the 
RNG  based  calculation  is  extended  to  include  the  O(b)  terms,  the  resulting  model  would 
perform  quite  well  in  comparison  to  other  second- order  closures. 

In  conclusion,  with  the  revised  coefficients  proposed  by  Yakhot  and  Smith3,  the  RNG 
K  —  e  model  now  performs  well  in  homogeneous  shear  flow  -  particularly  when  the  RNG 
based  anisotropic  eddy  viscosity  of  Rubinstein  and  Barton10  is  used.  The  RNG  second-order 
closure  model  needs  further  development,  however.  It  would  appear  that  an  extension  of 
the  rapid  pressure-strain  correlation  to  include  terms  of  0(b)  would  resolve  the  remaining 
deficiency  in  this  model.  Consequently,  our  current  assessment  of  RNG  based  turbulence 
models  is  now  more  optimistic  than  reported  earlier. 
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New  RNG 
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0 

-0.185 

0 

4.38 

Old  RNG 

K  —  e  Model 

0 

-0.489 

0 

11.70 

Standard 

K  —  e  Model 

0 

-0.217 

0 

4.82 

New  RNG 
Second- Order 
Closure 

0.489 

-0.091 

-0.244 

8.94 

Old  RNG 
Second- Order 
Closure 

0.533 

0 

-0.267 

00 

LRR  Model 

0.193 

-0.185 

-0.096 

5.65 

Experimented 

Data 

0.21 

-0.16 

-0.13 

4.8 

Table  1.  Comparison  of  the  equilibrium  values  of  the  various  models  with  the  experimental 
data  of  Tavoularis  and  Karnik7  on  homogeneous  shear  flow. 
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